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ABSTRACT 

We have used the full radial extent of images from the Hubble Space Telescope^s Advanced Camera 
for Surveys and Wide Field Planetary Camera 2 to extract surface brightness profiles from a sample 
of six, local lenticular galaxy candidates. We have modelled these profiles using a core-Scrsic bulge 
plus an exponential disk model. Our lenticular disk galaxies with bulge magnitudes My ^ —21.30 
mag have central stellar deficits, suggesting that these bulges may have formed from 'dry' merger 
events involving supermassive black holes while their surrounding disk was subsequently built up, 
perhaps via cold gas accretion scenarios. The central stellar mass deficits Mdof are roughly 0.5 to 2 
Mbh (black hole mass), rather than ^10 to 20 Mbh as claimed from some past studies, which is in 
accord with core-Sersic model mass deficit measurements in elliptical galaxies. Furthermore, these 
bulges have Sersic indices n 3, half light radii Re < 2 kpc and masses > IO^^Mq, and therefore 
appear to be descendants of the compact galaxies reported at z ^ 1.5 to 2. Past studies which have 
searched for these local counterparts by using single-component galaxy models to provide the z ^ 
size comparisons have over-looked these dense, compact and massive bulges in today's early-type disk 
galaxies. This evolutionary scenario not only accounts for what are today generally old bulges — which 
must be present in z ~ 1.5 images — residing in what are generally young disks, but it eliminates the 
uncomfortable suggestion of a factor of 3 to 5 growth in size for the compact, z 1.5 galaxies that 
are known to possess infant disks. 

Subject headings: galaxies: elliptical and lenticular, cD - galaxies: fundamental parameter - galaxies: 
nuclei - galaxies: photometry - galaxies: structure 
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1. INTRODUCTION 

It is now widely believed that all massive elliptical 
galaxies, and massive bulges in disk galaxies, house a 
supermassive black hole (SMBH) at their center (Magor- 
rian et al. 1998; Richstone et al. 1998; Ferrarese & Ford 
2005). In the standard cosmological formation paradigm 
of the universe, galaxies grow hierarchically, such that 
smaller systems merge to build larger ones (e.g.. White 
& Rees 1978; Khochfar & Burkert 2001). As galaxies 
containing SMBHs collide, their black holes will migrate 
to the center of the merger remnant through dynami- 
cal friction, and form a bound binary. The subsequent 
hardening of the black hole binary, following the dissi- 
pation of orbital energy and angular momentum to the 
nearby stars, impacts on the central stellar distribution 
of the newly merged galaxy. The gravitational slingshot 
ejection of these stars (from the binary's loss cone) via 
three-body interactions involving the binary is thought 
to be responsible for the physical origin of the partially 
depleted nuclear regions in luminous "core-Sersic" galax- 
ies which experienced 'dry', i.e. gas poor, merger events 
(Begelman et al. 1980; Ebisuzaki et al. 1991; Faber et al. 
1997; Milosavljevic & Merritt 2001; Merritt 2006; Sesana 
2010; Gualandris & Merritt 2012). 

Given the above scenario, the sizes and mass deficits of 
partially depleted cores are thought to reflect the amount 
of galactic merging and the ensuing extent of damage 
caused by binary SMBHs (after having eroded any preex- 
isting nuclear star clusters: Bekki & Graham 2010). Ac- 
curate measurements of the central stellar mass deficits 



in "core-Sersic" galaxies can therefore provide useful con- 
straints on models of galaxy formation and evolution. A 
key concept is that such core formation is by and large 
a cumulative process in which the ejected mass from 
SMBH binaries scales both with the final SMBH mass 
and the number of mergers (e.g., Milosavljevic & Merritt 
2001; Merritt 2006). Some recent studies have addition- 
ally postulated enhanced core depletions as a result of 
recursive core passages of recoiled SMBHs (Gualandris & 
Merritt 2008), or due to the actions of muhiple SMBHs 
from merging galaxies (Kulkarni & Loeb 2011). 

To date, several quantitative studies of the central stel- 
lar mass deficit have focused on elliptical galaxies (e.g., 
Graham 2004; Trujillo et al. 2004; Merritt 2006; Kor- 
mendy & Bender 2009; Dullo & Graham 2012). Vio- 
lent, major dry merger events are commonly, or at least 
traditionally, thought to produce these elliptical galaxies 
(Toomre & Toomre 1972; Barnes & Herquist 1992; KaufF- 
mann & Haehnelt 2000). Furthermore, as noted above, 
dry merger events involving supermassive black holes are 
also typically thought to be responsible for producing el- 
liptical galaxies with partially depleted cores. Therefore, 
some ambiguity exists regarding the mechanisms for the 
formation of the lenticular (SO) disk galaxies having de- 
pleted cores. 

In 1936 Hubble introduced lenticular galaxies into the 
tuning fork diagram (Jeans 1928; Hubble 1929); they 
were added as a hypothetical transition type between 
elliptical and spiral galaxies. Later observation of the 
populations of lenticular and spiral galaxies in clusters re- 
vealed the dominance of lenticular galaxies (spiral galax- 
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ies) in local (in distant) clusters (e.g., Dressier 1980; 
Dressier ct al. 1997; Couch et al. 1998; Fasano at al. 
2000; Desai ct al. 2007; Poggianti et al. 2009; Sil'chenko 
et al. 2010). This suggested an evolutionary transfor- 
mation of SO galaxies from spiral galaxies, where mecha- 
nisms such as ram pressure stripping (Gunn & Gott 1972; 
Quilis et al. 2000), galaxy- galaxy harassment (Moore et 
al. 1996), strangulation (Larson et al. 1980; Bekki et al. 
2002; Kawata & Mulchaey 2008) and gas starvation by 
AGN (van den bergh 2009) were forwarded to account 
for the removal of disk gas and the subsequent quench- 
ing of star formation, thereby passively fading the spiral 
galaxies and erasing their spiral arms. These formation 
scenarios, however, do not in themselves explain the cen- 
tral stellar depletions observed in some SO galaxies. 

The theoretical work of Steinmetz & Navarro (2002) 
has suggested that a galaxy's morphology is a transient 
phenomenon. In this hierarchical picture, classic ellip- 
tical galaxies are built through major mergers (of disk 
galaxies) and may progressively regrow stellar disks, by 
gas accretion, which remain intact only until the next sig- 
nificant merger (see Okamoto & Nagashima 2001; Gov- 
ernato et al. 2009; Pichon et al. 2011; Sales et al. 2011; 
Conselice et al. 2012 for supporting arguments). The 
number of such cycles may however be low (i.e. 1) rather 
than several (3 to 5). Likewise, Arnold et al. (2011) and 
Forbes et al. (2011) recently pointed out that SOs might 
form through a two-phase inside-out assembly with the 
inner regions built early via a violent major merger, and 
"wet" minor mergers subsequently contributing to the 
outer parts. 

This assembly scenario is consistent with the observed 
presence of partially depleted cores in luminous SO galax- 
ies. It also suggests that their bulges, most of which we 
know are old (e.g., MacArthur, Gonzalez & Courteau 
2009, and references therein), were already around at 
z K, 1.5 to 2. Graham (2013, his Fig. 1) revealed that 
the compact galaxies at z = 1.5 to 2 have masses and 
structural properties consistent with those of the brighter 
bulges in local disk galaxies. Rather than being the pre- 
cursors of elliptical galaxies prior to significant size evolu- 
tion, Graham (2013) advocated that some of these com- 
pact high-z galaxies may therefore be associated with 
the bulges of modern disk galaxies. Furthermore, in this 
scenario in which compact galaxies evolve into disk galax- 
ies with compact bulges, the (high velocity)-end of the 
galaxy 'velocity dispersion function' would not be ex- 
pected to evolve from z ^ 1.5 to 0, just as observed 
(Bezanson, van Dokkum & Franx 2012). 

Regarding the analysis of surface brightness profiles, 
Trujillo et al. (2004) avoided the inclusion of SO galaxies 
in their sample of well-resolved local galaxies because it 
was felt that the introduction of an 8-parameter model 
(5 or 6 core-Scrsic parameters for the bulge plus 2 for the 
exponential disk) might be considered excessive at that 
time. Ferrarese et al. (2006) also avoided introducing 
two additional disk parameters in their modelling of well- 
resolved early- type galaxies in the Virgo cluster; however 
because their sample included disk galaxies, their core- 
Sersic parameters do not correspond to the bulge compo- 
nenlQ. Although Kormendy et al. (2009) did fit a Sersic 

^ The outer flattened disk of a lenticular galaxy, if not accounted 
for, can effectively result in a single 'galaxy' Sersic index which is 



bulge plus an exponential disk to their SO galaxies, they 
marked the core region by eye rather than objectively 
fitting a core-Sersic bulge plus exponential disk model. 
Depending on the sharpness of the (inner core)-to-(outer 
Sersic) transition region, this practice may substantially 
over-estimate the formal 'break radius' — which is not to 
be confused with the inner or outer edge of the transition 
region. 

Here we endeavor to provide the most accurate mea- 
surements to date of the centrally depleted stellar mass 
in lenticular disk galaxies. This is achieved by simulta- 
neously fitting both a core-Scrsic model to the bulge and 
an exponential model to the disk. 

This paper is organized as follows. Our initial sample 
of six suspected lenticular galaxies with depleted cores, 
plus the data reduction, and the light profile extraction 
technique are discussed in Section [2] In Dullo & Graham 
(2012) we did not have a sufficient radial extent of these 
six galaxies' light profiles to include a disk component. 
Here we are able to test if this may have slightly biased 
some of our previous measurements of the bulge param- 
eters. Section [3] introduces the analytical models used 
to describe the light profiles of "core-Sersic" lenticular 
galaxies, and Section |4] details our fitting analysis while 
in Section [5] we describe the method for measuring the 
mass deficits of the core-Sersic lenticular galaxies. Our 
interpretations of the central stellar mass deficits in the 
context of core-Sersic lenticular galaxy formation scenar- 
ios are discussed in Section [51 In Section 16.11 we discuss 
the formation of core-Scrsic lenticular galaxies, and in 
Section 16.21 we compare the physical properties of their 
bulges with those of compact galaxies at z = 1.5 to 2. 
In Section [^31 we discuss the role of galaxy environment. 
Section [7] summarizes our main conclusions. At the end 
of this paper are two appendices. The first provides notes 
on the six individual galaxies studied in this paper and 
the second provides a response to the criticism of Dullo 
& Graham (2012) by Lauer (2012) regarding the identi- 
fication of partially depleted cores. 

2. DATA AND PHOTOMETRY 

2.1. Sample selection 

We have targeted the six 'core-Sersic' lenticulaiQ galax- 
ies from the sample of 39 relatively bright, nearby, early- 
type galaxies analysed by Dullo & Graham (2012). This 
initial sample of 39 galaxies came from Lauer et al. (2005) 
who had claimed that they all have partially depleted 
cores. Using the core-Sersic model rather than the Nuker 
model, Dullo & Graham (2012) subsequently revealed 
that seven of them did not have partially depleted cores 
relative to the inward extrapolation of their outer Sersic 
profiles (see also l9.2p . The six, suspected lenticular galax- 
ies with depleted cores are listed in Table [TJ 

2.2. Imaging data 

High resolution Hubble Space Telescope (HST) opti- 
cal images of these galaxies were retrieved from the 
public HST data archive. We only used images from 

larger than that of the triaxial bulge component (e.g., Meert et 
al. 2012). Use of this larger 'galaxy' Sersic index can result in an 
inflated measurement of the bulge's central mass deficit. 

^ The initial morphological classification was taken from the 
Third Reference Catalog, FIC3, de Vaucouleurs et al. (1991). 
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broad-band filters. While the HST Advanced Camera 
for Survey (ACS; Ford et al. 1998) Wide Field Chan- 
nel (WFC) F475W and F850LP images, from the Virgo 
Cluster Survey GO-9401 program (PI: P. Cote), are 
available for NGC 4382, for the remaining five galax- 
ies (NGC 507, NGC 2300, NGC 3607, NGC 3706 and 
NGC 6849) we used the Wide Field Planetary Cam- 
era 2 (WFPC2; Holtzman et al. 1995) F555W images 
taken from the following programs: NGC 507, NGC 
3706 and NGC 6849 from the GO-6587 program (PI: 
D. Richstone); NGC 2300 from the GO-6099 program 
(PI: S. Faber) and NGC 3607 from the GO-5999 pro- 
gram (PI: A. Phillips). For NGC 3607, we also used the 
near-infrared NICM0S/NIC2 F160W image observed in 
the GO-11219 program (PI: A. Capetti) to enable us to 
better avoid this galaxy's inner dusty spiral structure. 
Global properties and observation details of all the sam- 
ple galaxies are listed in Table [TJ 

The WFPC2 consists of 3 wide field cameras (WF2, 
WF3 and WF4) and a high resolution Planetary Camera 
(PCI); each has a CCD detector with 800x800 pixels. 
The three wide field cameras have a 0".l per pixel spatial 
sampling. The smaller, high resolution planetary camera 
(PCI), where each galaxy's center had been placed, has 
a plate scale of 0".046 and a square 36"x36" field of view 
(FOV). 

Having two CCDs cameras with 2048 x 4096 pixels, 
the ACS Wide Field Channel has a 0".049 pixel scale 
and covers a 202" x 202" rhomboidal area. 

Combining the constituent CCDs of each camera, the 
light profiles from the WFPC2 and ACS observations 
probe a large range in radius, R > 80". For comparison, 
Lauer et al. (2005, and references therein) deconvolved 
all six sample galaxy images with the PSF and provided 
profiles out to a maximum of 10". The greater radial ex- 
tent that we now have, from high quality CCD imaging, 
enables us to better sample and measure the disk light. 

While our WFPC2/F555W (roughly Jonson-Cousins 
y-band) images match Lauer et al.'s. (2005) F-band 
images, for NGC 4382 the ACS/F475W (roughly SDSS 
g) and ACS/F850LP (roughly SDSS z) data arc trans- 
formed to the V-band using 

g-y = 0.98 (g-2)-1.43, (1) 

which is acquired here from the least squares fit to the 
{g — V) and (g — z) data shown in Fig. [1] see also Sirianni 
et al. (2005) and Kormcndy et al. (2009). The zero points 
in the Vega magnitude systems, which are adopted here 
to calibrate the WFPC2 (Hohzman et al. 1995) and ACS 
(Sirianni et al. 2005) profiles, were taken from the STScI 
web sit(j3- 



2.3. Data Reduction 

The WFPC2, ACS and NICMOS images, retrieved 
from the Hubble Legacy ArchivcQ (HLA), were processed 
using the standard HLA reduction pipeline. The reduc- 
tion steps include bias subtraction, geometric distortion 

^ The ACS/WFC images used here are taken from observations 
before 2006. 

* http://hla.stsci.edu 




1.95 2 2.05 

g - z 

Fig. 1. — Photometric transformation from ACS g and 2-band 
magnitude systems into WFPC2 l/-band magnitude for NGC 4382. 
We use our g and 2:-band plus Lauer et al.'s (2005) V-band sur- 
face brightnesses of the galaxy at different radii to obtain the data 
points. The line is the least-squares fit to the data (Eqn.[l}. 



TABLE 1 

Global properties for our sample of six suspected 'core-Sersic' 
lenticular galaxies. 



Galaxy 


Type 




D 


C7 


Filter 


Exp. Time 






(mag) 


(Mpc) 


(km s-i) 




(s) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


NGC 0507 


SO 


-22.56 


63.7" 


306 


F555W 


1700 


NGC 2300 


SO 


-21.33 


25.7" 


261 


F555W 


1520 


NGC 3607 


SO 


-21.55 


22.2' 


224 


F555W 


160 












F160W 


1151 


NGC 3706 


SO 


-22.08 


45.2" 


270 


F555W 


1400 


NGC 4382 


SO 


-21.38 


17.9' 


179 


F475W 


750 












F850LP 


1210 


NGC 6849 


SBO 


-22.51 


80.5" 


209 


F555W 


900 



Notes. — Col. (1) Galaxy name. Col. (2) Morphological type from RC3 
(de Vaucoulers et al. 1991). Col. (3) Absolute V-band bulge magnitude 
(galaxy magnitude for NGC 3706 since wc (ultimately) adopt an elliptical 
morphology, see Section 4.1) obtained using our bulge-to-disk (B/D) flux 
ratios and absolute y-band galaxy magnitudes from Lauer et al. (2007b). 
These magnitudes are corrected for Galactic extinction, inclination and 
internal dust attenuation (Driver et al. 2008, their Table 1 and Eqs. 1 
and 2) and (1 -|- z)* surface brightness dimming, and adjusted using the 
distance from col. (4). Sources: (t) Tonry et al. (2001) after reducing 
their distance moduli by 0.06 mag (Blakeslee et al. 2002); (n) from 
NED (3K CMB). Col. (5) Central velocity dispersion from HyperLcdifl 
(Paturel et al. 2003). 
''l |http://leda.univ-lyonl.fr| 

correction, dark current subtraction and flat fielding us- 
ing date-appropriate references: the WFPC2 (McMaster 
et al. 2008) and ACS (Pavlovsky et al. 2006) data and 
instrument handbooks provide detailed descriptions of 
these steps. 

The automatic HLA reduction pipeline also subtracts 
the sky background values from the images. Contamina- 
tion from a poor subtraction is a possibility which makes 
the automatic sky background estimation of the spatially 
extended galaxies slightly uncertain, it essentially affects 
the outermost parts (fainter regions) of such galaxies' 
profiles. Having major-axis diameters of < 4' and minor- 
axis diameters of < 3' (using fis^'^^ mag arcsec"^ as 
a reference level, NED), NGC 507, NGC 2300, NGC 
3706 and NGC 6849 are within the WFPC2 field of view, 
thus, they are less prone to the sky subtraction errors. 
NGC 3607 and NGC 4382, however, extend beyond the 
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Fig. 2. — Major-axis surface brightness, elliptieity, position angle (measured in degrees from north to east) and isophote shape parameter 
(B4) profiles for the galaxies in Table [T] Top right panels show the surface brightness profiles obtained with the F555W (~ V-band) filter, 
with the exception that a NICMOS NIC2 F160W (i?-band) profile is also shown for NGC 3607, and NGC 4382 was imaged in the ACS 
F475W (g-band) and F850LP (z-band) filters. 
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WFPC2 and ACS field of views, respectively. Ferrarese 
et al. (2006) used a combination of the ACS and ground- 
based data to better constrain the sky level around NGC 
4382, this data was later re-modelled by DuUo & Graham 
(2012, their Fig. 7). We find here that our light profile, 
and fit parameters for this galaxy, match those in DuUo 
& Graham (2012), indicating a small or zero sky subtrac- 
tion error by the pipeline. This gives us some confidence 
that the outermost profile of NGC 3607 is also correct. 

Finally, the images were masked to avoid the gaps be- 
tween individual CCD detectors, the partially missing 
quadrant of the WFPC2 images, and regions including 
bright foreground stars, background galaxies, chip de- 
fects and galaxy dust lanes. 

2.4. Surface photometry 

The galaxy light profiles were extracted using the 
IRAF/STSDAS task ellipse (Jcdrzejewski 1987). We 
used the ellipse task to construct the best-fitting con- 
centric elliptical isophotes to the galaxy image starting 
from an initial elliptical isophote defined by the first 
guess values of the isophote center (X,Y) , ellipticity (e) 
and position angle (PA). We used median filtering, log- 
arithmic spacing, and 3a clipping for flagging deviant 
sample points at each isophote. To extract the major- 
axis light profiles, for all galaxies, the isophote centre (X 
,Y), e and PA were set free to vary. For each galaxy we 
found that the isophote center was stable, i.e. within the 
small quoted error box. 

The deviations of isophotes from perfect ellipses can 
be characterized well by higher order coefficients of the 
Fourier series expansion of the intensity (Jcdrzejewski 
1987). In particular, the coefficient of the cos 49 term, 
i?4, is found to be of great importance in describing these 
deviations. Negative B4 values imply that the isophotes 
are "boxy", but, if B4 is positive, the isophote will be 
identified as "disky" . Disky isophotes are frequently due 
to the presence of embedded disks in less massive, fast ro- 
tating, dissipative (gas rich) galaxies (e.g.. Carter 1978, 
1987; Davies et al. 1983; Bender et al. 1988; Peleticr 
et al. 1990; Jaffe et al. 1994; Faber et al. 1997). In 
contrast, dissipationless violent relaxation of stars are 
often invoked to explain the boxy isophotes of triax- 
ial, pressure-supported massive galaxies (Nieto & Bender 
1989), which usually also contain hot X-ray emitting gas 
(e.g., Cattaneo et al. 2009). 

In Fig. [2] we show the results of the ellipse fitting. 
As noted in Rest et al. (2001), the isophotal parameters 
derived from the ellipse fits tend to be more uncer- 
tain in the very inner regions (i? < 0."3); this feature 
is mainly due to the PSF and some contributions from 
the discrete sampling and subpixel interpolation of the 
IRAF fitting routine (Ravindranath et al. 2001). We 
note that while the disky/boxy isophotal deviations vary 
with radius in all of the six galaxies, NGC 507 and NGC 
4382 have core regions which appear preferentially boxy 
within 1" (Fig.H]). Nieto & Bender (1989) reported that 
the SO galaxy NGC 2300 has boxy inner (1" < R< 20") 
isophotes and disky 'pointed' outer (20" < R < 50") 
isophotes, in agreement with our result (Fig. 2). 

In constructing 'composite' light profiles for each 
galaxy, we combine the very inner (i? < 1") Lauer et 
al. (2005) PSF-deconvolved, high resolution y-band light 
profile with our photometrically calibrated (1^-band) pro- 



file of the PSF-unaffected region beyond ^ 1". For all 
the six galaxies in our sample, we find an excellent agree- 
ment (an overlap) between our light profiles and those 
published by Lauer et al. (2005) over the 1" — 5" radial 
range. 

3. MODELS FOR 'CORE-SERSIC LENTICULAR GALAXIES 

Given the two component (bulge/disk) nature of lentic- 
ular galaxies, we apply a bulge-to-disk photometric de- 
composition to the one-dimcnsional, major-axis surface 
brightness profiles of all the galaxies in our sample. 
When needed, a bar component is also included. 

The radial intensity distribution of our lenticular 
galaxies' disk component is modelled with an exponential 
function given by 



Idisk{R) = hd^^V [-R/h], 



(2) 



where /od and h are the central intensity and scale length 
of the disk respectively. 

We adopt the Ferrers (1877) function to describe the 
radial intensity distribution of the bar component, given 
by 

Ibar{R) = lobar [l — {R/o,bar)^] , ^ < a&ar, 

har{R) =0, R> ahar, 

(3) 

where lobar, i&ar and n\jar srs the central intensity, 
major-axis length and shape parameter of the bar, re- 
spectively (cf. Laurikainen et al. 2010). 

Since the work by Caon et al. (1993) and Andredakis, 
Peleticr & Balcells (1995), several studies revealed that, 
in general, the Sersic (1963) model describes the under- 
lying light distributions of both elliptical galaxies and 
the bulges of disk galaxies exceedingly well over a large 
radial range. This model can be written as 



I{R) — lo exp 



(4) 



where Iq = /(i? = 0) is the central intensity. The quan- 
tity bn « 2n- 1/3, for 1 < n < 10 (e.g., Caon et al. 1993; 
Graham 2001) is a function of the Sersic index n, and is 
defined in such a way to ensure that the half light ra- 
dius. Re, encloses half of the total luminosity. Reviewed 
in Graham & Driver (2005), the total luminosity of the 
Sersic model within any radius R is given by 



LT,Seri< R) = leRpirn 



(^7(2n,x), (5) 



where ^(2n,x) is the incomplete gamma function and 

X = bn{R/Rey/". 

Systematic downward departures of the inner light pro- 
file relative to the inward extrapolation of the outer 
Sersic model are known to exist in luminous galaxy /bulge 
light profiles. These are not exactly the same objects as 
"core" galaxies identified with the Nuker model (Lauer 
et al. 1995; Byun et al. 1996) or the double power-law 
model of Ferrarese et al. (1994). "Core-Sersic" galaxies 
have partially-depleted cores relative to their outer Sersic 
profile whereas "core" galaxies need not have any deficit 
but instead simply an inner power-law slope 7 < 0.3. 
As detailed in DuUo & Graham (2012), the core-Scrsic 
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model (Graham et aL 2003), a combination of an inner 
power-law and an outer Sersic function, provides a good 
representation of the brightness profiles of bulges in disk 
galaxies with depleted cores. The model is defined as 



I{R) = /' 
with 



7/c 



exp 



-b„ 



Q I pa \ l/{an) 



R°' +R 



(6) 



/' = /b2-^/" exp 6„(2i/"i?fc/i?e) 

(7) 

The term /{, denotes the intensity at the core's break 
radius i?h, 7 is the slope of the inner power law, and a 
controls the sharpness of the transition between the inner 
power-law and the outer Sersic profile. Both Rg and b„ 
have the same general meaning as in the Sersic model. 

4. FITTING ANALYSIS 

Fig. [3] displays the best model fits to the major-axis 
surface brightness profiles of the 6 galaxies listed in Ta- 
ble [TJ We follow the iterative fitting technique of Dullo & 
Graham (2012), minimizing the root mean square (rms) 
residuals to determine the best fit parameters that match 
the data. It is important to note that these profiles 
now cover a large range in radius (out to i? > 80"), 
giving both the core-Sersic and the exponential models 
enough radial expanse to both quantify the curvature 
in the bulge profile (i.e. the Sersic index) and define 
the disk scale length. Table [2] lists the fit parameters 
obtained from the adopted core-Sersic plus exponential 
bulge-l-disk models. 

Apparent from Fig. [3] is that the global light distri- 
bution of all the sample galaxies, except for NGC 3706 
(E not SO) and NGC 6849 (barred), can be accurately 
represented by the core-Sersic -|-exponential model with 
a small rms residuals < 0.05 mag arcsec"^. If a single 
core-Sersic model is fit to the light profile of a disk galaxy 
(i.e. a two-component system comprised of a flattened 
disk and a triaxial bulge), this would result in parame- 
ters that do not describe the triaxial bulge component of 
the galaxy, effecting in particular the Sersic index, and 
ultimately the mass deficit measurement as happened in 
Ferrarese et al. (2006, See also Meert et al. 2012). Al- 
though NGC 3706 is classified as a lenticular galaxy in 
the RGB, we find that its light profile is best fitted by 
a single core-Sersic model, without an exponential disk 
component. This suggests that the galaxy may be an 
elliptical galaxy misclassified as an SO in the RG3, con- 
sistent with the conclusion of Laurikainen et al. (2010). 
As for NGG 6849, the 3-componcnt bulge-bar-disk light 
profile is well described by the core-Sersic bulge-l-Fcrrers 
bar+exponential disk decomposition model (see Fig|3] 
and Section r9.1.6|) . 

While the F-band surface brightness profile of 
NGC 3607 is well described by the core-Sersic 
bulge-|-exponential disk model, the 1".31 core radius of 
this fit is suspiciously larger than the 0".22 core radius 
reported by Richings et al. (2011) from their i/-band 
profile analysis. The explanation is that a dusty nuclear 
spiral has caused a reduction to the inner y-band light 
profile. Our analysis of the PSF-affected i7-band light 
profile (Fig. |3]) gives a core radius of 0".ll, which is in 



fair agreement with the result found by Richings et al. 
(2011) but does not support a large core in this galaxy. 
Given that the core size we find with the PSF-affected 
iJ-band data is comparable to the seeing, coupled with 
the dusty spiral structure, we arc unable to conclude if 
there is indeed any partially-depleted core in this galaxy, 
and as such we exclude it from our final analysis. 

We have found that the bulge model parameters (Ta- 
ble [5]) generally agree well with those from Dullo & Gra- 
ham (2012). For NGC 2300, the contribution of the disk 
light to the inner 15" profile which we modelled in Dullo 
& Graham (2012) did however result in a break radius 
Rb = 0".98 which is slightly larger than the one from this 
work {Rh ~ 0".53). With the exception of this galaxy, 
the agreement between the break radii from these two 
studies is good, constrained to less than a 20% discrep- 
ancy. For each galaxy, we find that the y-band surface 
brightness at the break radius (yUb.y) given in Table 2 
agrees with those presented in Dullo & Graham (2012) 
to within the range of error bar: there is ~ 5% discrep- 
ancy, which is well within the Dullo & Graham (2012) 
^ 10% uncertainty range. We have also compared the in- 
ner power- law slopes (7) from the two studies. We find 
that the new 7 values (Table 2) are consistent with those 
reported in Dullo & Graham (2012). 

The five core-Sersic lenticular galaxies (Table [5]) have 
bulge Sersic indices n ^ 3, while typically for core-Sersic 
elliptical galaxies n is greater than 3. This may be be- 
cause, for a given magnitude, the bulges of disk galaxies 
appear to have lower Sersic indices than elliptical galax- 
ies of the same magnitude (Graham 2001, his Fig. 14). It 
may also be that the elliptical core-Sersic galaxies from 
Dullo k Graham (2012) are brighter. 

5. CENTRAL STELLAR MASS DEFICIT 

The central stellar mass deficit (Mdof) of 'core-Sersic' 
galaxies is predicted to be generated mainly via the three- 
body encounters of the core stars with the inspiraling 
black holes in a merger remnant (e.g., Ebisuzaki et al. 
1991; Milosavljevic & Merritt 2001; Merritt 2006). Gra- 
ham (2004) measured this stellar mass deficit from the 
luminosity difference idef between the inward extrapola- 
tion of the outer Sersic model and the (sharp-transitioij^ 
core-Sersic model. This approach was adopted in subse- 
quent works by Ferrarese et al. (2006), Merritt (2006) 
and Hyde et al. (2008). We adopt a slightly different 
methodology to those studies by using a finite value for 
a in Eq. [6] and thus a smoother transition region be- 
tween the inner power-law and the outer Sersic profile 
of the core-Sersic model. The total core-Sersic model 
luminosity (Trujillo et al. 2004; their Eq. A19) is given 
by 

Lt.,s = 2^I'n{Re/b^f J 



X a; 



[b'^Rb/ReY 



■i(2-7-a)/a 



dx. 



where all the parameters have the same meaning as in 
Eq. [6l The difference in luminosity between the outer 

^ Trujillo ct al. (2004) set the transition parameter « — > 00 to 
obtain the 5-parameter sharp-transition core-Sersic function given 
by their Eqn. 5. 
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TABLE 2 

Structural parameters. 



Galaxy 


Type 




Rb 


Rb 


7 


a 


n 


Re 


Re 




h 


(B/T)obs 


(B/T)cor 








(arcsec) 


(pc) 








(arcsec) 


kpc 




(arcsec) 






n 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


NGC 0507 


SO 


16.38 


0.33 


102 


0.07 


5 


2.19 


5.34 


1.65 


21.03 


27.69 


0.32 


0.43 





NGC 2300 


so 


16.61 


0.53 


70 


0.08 


2 


2.20 


7.69 


1.02 


20.39 


21.08 


0.42 


0.57 


44 


NGC 3607 


so 












2.39 


7.75 


0.84 


19.16 


18.61 


0.40 


0.57 


59 


NGC 3706t 


E 


14.16 


0.11 


24 


-0.02 


10 


6.36 


42.08 


9.18 










53 


NGC 4382 


SO 


15.01 


0.27 


24 


0.07 


5 


2.65 


11.14 


0.99 


19.50 


35.07 


0.28 


0.40 


39 


NGC 6849* 


SBO 


16.67 


0.18 


69 


0.20 


5 


3.23 


7.77 


2.98 


20.72 


16.93 


0.31 


0.46 


55 



Notes. — Structural parameters from fits to the y-band major-axis surface brightness profiles (Fig. [S]!. Col. (1) Galaxy name. Col. (2) 
Our adopted morphological classification. Col. (3)-(10) Best-fit bulge structural parameters from the core-Sersic model, Eq. [B] Col. 
(11)-(12) Best-fit disk structural parameters from the exponential model, Eq. |3] Break surface brightness lJ,b,v a-iid disk central surface 
brightness /xod v are in mag arcsec"^. Col. (13) Expected V-band bulge-to-total (B/T)obs flux ratios obtained using Graham & Driver 
(2005, their Eq. 19). These ratios are corrected for Galactic extinction, surface brightness dimming, internal dust attenuation and 
inclination (Driver et al. 2008, their Table 1 and Eqs. 1 and 2) and listed in Col. (14). Col. (15) Disk inclination angles (i) arc derived 
using the galaxies' major- and minor-axis diameters from NED. f We classify NGC 3706 as an elliptical galaxy based on our light profile 
analysis. | The fit parameters for the bar component: /xoi,ar,v=21.63 mag arcsec"^, ai,ar = 29". 04, and ni,ar = 6.37. 
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Fig. 3. — Fits to the major-axis surface brightness profiles of the galaxies in Table [Tl For NGC 3607 (with a dusty nuclear spiral), the 
core-Sersic model fit to the PSF-affectod _ff-band data (crosses) is shown in addition to the core-Sersic-|-exponential fit to the F-band data 
(dots). The dashed curves show the Sersic component of the core-Sersic fits to the data, while the dotted curves show the exponential 
function simultaneously fit to the large scale disks. The dash-dot-dot-dot curve shows the Ferrers function fit to the bar component of 
NGC 6849. The solid curves show the complete fit to the profiles, with the rms residuals. A, about each fit given in the lower panels. For 
NGC 3607, we only show the rms residual from the K-band profile fit, while A from its _f/-band profile fit is 0.014 mag arcsec^^. Data points 
excluded from the fits arc shown by the open circles (see Section |9ll. 
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Fig. 4. — Central mass deficit {M^^{) versus black hole mass 
(Mgu) for 'corc-Sersic' galaxies listed in Table [2] and |3] The only 
elliptical galaxy is circled. The M-a relation presented in Gra- 
ham ct al. (2011) was used for estimating the SMBH masses of 
the galaxies. A horizontal solid line connects the two predicted 
black hole masses of NGC 6849: i) assuming it is a barred galaxy 
(triangle) and ii) considering it as a non-barred galaxy (star). A 
representative error bar is shown at the bottom of the panel. 

Sersic model (Eq. [S]) and the core-Scrsic model (Eq. [S]) is 
of course the central stellar flux deficit. 

To convert the luminosity deficits into mass deficits, in- 
dividual stellar mass-to-light [M/L) ratios were derived 
for each galaxy. To determine these ratios, we have made 
use of the available (preferentially nuclear) colours of the 
galaxies as well as the colour-age- metallicity- (M/L) dia- 
gram from Graham & Spitler (2009; their Fig Al). As in 
Graham (2004), we have assumed that the inner regions 
of these galaxies have an evolved (old) single population 
of stars. 

For NGC 2300, NGC 3607 and NGC 4382, the nuclear 
V — I colours of 1.33, 1.36 and 1.11, taken from Lauer 
et al. (2005), correspond to stellar M/Ly ratios of 5.0, 
5.7 and 2.6, respectively. Using the HyperLeda database, 
NGC 507 and NGC 3706 have V -I colours of 1.40 and 
1.34, respectively, while NGC 6849 has ^-7=1. 03. From 
this, we obtain M/Ly ratios of 5.5, 5.2, and 2.4 for NGC 
507, NGC 3706 and NGC 6849, respectively (see Table 
[3]). These mass-to- light ratios were used to convert the 
stellar flux deficits into stellar mass deficits. These have 
been plotted in Fig. U] against each galaxy's expected 
black hole mass as derived from the M — a relation (Gra- 
ham et al. 2011; the lower half of their Table 2) based on 
the a values in Table [TJ 

As noted above, high-accuracy N-body simulations by 
Merritt (2006) revealed that the effect of multiple dissi- 
pationless mergers on core formation is cumulative, thus, 
the extent of core evacuation could reflect the amount 
of merging. He found that the total mass deficit after 
N 'dry' major mergers is « 0.57VMbh, where Mbh is 
the final black hole mass. This result hinges on a key 
assumption that the central black hole would tidally dis- 
rupt infalling high-density satellites, thereby protecting 
any preexisting core (Fabcr ct al. 1997; Boylan-Kolchin 
& Ma 2007). 

Fig. 0] shows that the mass deficits for the 'corc-Scrsic' 
lenticular galaxies (plus one elliptical galaxy) listed in 



Table [5] are Mdef ~ 0.5 — 2A/bh, which translates to a 
few (1 to 4) 'dry' major merger events. For reference, 
measurements of close galaxy pairs in real data (i.e. not 
simulations) suggest that massive galaxies (> IO^^'^Mq) 
have experienced 0.5 to 2 major mergers since z 0.7 
to 1 (e.g., BeU et al. 2004, 2006; Bluck ct al. 2012; Man 
et al. 2012; Xu et al. 2012). Our figure is double this 
value, but it has no upper redshift constraint associated 
with it. The mean elliptical galaxy Mdcf/MeH mass ra- 
tio from Graham (2004) is 2.1 ±1.1, while the mean ratio 
from Ferrarese et al. (2006) is 2.4 ± 0.8 after excluding 
the SO galaxy NGC 4382 from their sample (as they did 
at the end of their Section 5.2). We used the Mbh-o" re- 
lation presented in Graham et al. (2011; the lower half of 
their Table 2) for estimating the SMBH masses of the five 
core-Sersic galaxies shown in Table [31 While this Mbh- 
(T relation was constructed by combining corc-Scrsic and 
Sersic galaxies, Graham & Scott (2012, their Table 3) 
have recently shown that core-Sersic and Sersic galax- 
ies follow similar AfsH-c relations. They additionally re- 
ported an AfBH-f relation for their combined (core-Sersic 
+ Sersic) galaxies which is consistent with the relation 
found by Graham et al. (2011). We find that, for all our 
core-Sersic galaxies, the SMBH masses predicted using 
the Graham & Scott (2012) M-Q-^-a relations which are 
defined by the core-Sersic, Sersic or core-Scrsic + Sersic 
galaxies agree with those used in this study (Table 3) 
within the range of error bars. We make use of Eq. 4 
from Graham et al. (2011) and the cr-values in Table [1] 
and assume a 10% uncertainty on a to estimate the la 
uncertainty on each galaxy's SMBH mass (Table [3]). 

For NGC 4382, Kormendy & Bender (2009) reported 
a mass deficit of Mdcf ^ 1-3 x 1O^M0 w I'UMbh when 
using Mbh = 1 x IO^A/q, while Giihckin et al. (2011) 
used Nuker model parameters and found a mass deficit 
of Afdcf ^ 5.9 X lO^A^Q « 45.6A//bh when they assumed 
a very small black hole mass of 1.3 x IO^ALq. Ferrarese 
et al. (2006, see their Fig. 20) noted that the large-scale 
stellar disk of NGC 4382, which was incorporated into 
their single Sersic fit to this galaxy, might have biased 
their ~ \QMbh mass deficit measurement. Indeed, their 
galaxy Sersic index was « 6-7, while we have found that 
the bulge Sersic index is less than 3. Our analysis of this 
bulge's central stellar mass deficit yields Afdef ^ 1-2 x 
IO^ALq « 1.3A/BH (using AfsH = 9.55 x lO'^Afg), smaller 
than the results of Ferrarese ct al. (2006), Kormendy et 
al. (2009) and Giitekin et al. (2011). 

Previously, Milosavljevic & Merritt (2001), Milosavl- 
jevic et al. (2002), Ravindranath et al. (2002) and Lauer 
et al. (2007a) had used the Nuker model to measure mass 
deficits that are up to an order of magnitude larger than 
those obtained here. As detailed in Trujillo et al. (2004) 
and Dullo & Graham (2012), this discrepancy arises in 
part from the differences in the estimated core sizes from 
the two models. While the Nuker model break radii are 
up to 3 times larger than the core-Sersic model break 
radii, due to the Nuker model fitting a power-law to each 
galaxy's outer curved Sersic profile (Graham ct al. 2003), 
the core-Sersic break radii agree with model-independent 
core sizes where the negative logarithmic slope of the 
fight profile equals 0.5 (Dullo & Graham 2012). 

More recently, Kormendy et al. (2009) fit Sersic models 
to core-Sersic galaxy light profiles and tried to quantify 
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the core from a visual inspection, rather than using the 
core-Sesic model in an objective analysis. They report 
large break radii and mass deficits Mdof 10 — 2QMbh- 
In an effort to reconcile these large deficits, they specu- 
lated that they may be due to recoiled black holes, which 
might enhance the core depletion following their repeti- 
tive core passages (Boylan-Kolchin et al. 2004; Gaulan- 
dris & Merritt 2008); cf. also Postman et al. (2012) for 
a similar reasoning. If this mechanism always occurred, 
then the lower Mdci/MsH ratios in Dullo & Graham 
(2012) for elliptical galaxies would suggest that the effec- 
tive number of major mergers is less than one, i.e. only a 
minor merger is required. However given the near-linear 
relation between black hole mass and host spheroid mass, 
and luminosity, for core-Sersic galaxies (Graham 2012; 
Graham & Scott 2012) this scenario is unlikely. At least 
for the elliptical core-Scrsic galaxies, which are the dom- 
inant population among the known core-Scrsic galaxies, 
this linear black hole-galaxy mass scaling relation readily 
arises from the self addition of comparable mass systems, 
not minor mergers. 



6. DISCUSSION 

6.1. Formation of core-Sersic lenticular galaxies 

As noted above, the central stellar mass deficits for the 
'core-Sersic' lenticular galaxies are found to be Mdef ~ 
0.5 — 2Mbh (FigH]). By comparing these mass deficits 
with the results of Merritt (2006), one would conclude 
that the bulges of core-Scrsic SOs have experienced the 
equivalent of a few 'dry' major merger events. Devia- 
tions from this rule may however exist. For example, 
substantial black hole recoiling events would lower this 
figure of a few, while 'loss cone' re-filling through wet 
mergers and the production of new stars would allow 
scope for an increased number of mergers. This latter 
scenario is unlikely given, in general, the red colours 
of the galaxies' bulges. However, Carter et al. (2011) 
have found FUV — NUV colour gradients in several 
core-Scrsic galaxies which seem to be inconsistent with 
the major, dry merger scenario detailed in Faber et al. 
(1997). Also, the globular cluster specific frequencies of 
some intermediate luminosity ellipticals are found to be 
systematically lower than the ones in bright, core-Sersic 
elliptical galaxies (Harris & van den Bergh 1981; van den 
Bcrgh 1982). This may reflect that core galaxy formation 
may not involve dry merger events, because such merg- 
ers are thought to be inefficient at creating new globular 
clusters. 

Although not (yet) popular, various studies have pro- 
posed alternative mechanisms for generating cores in lu- 
minous galaxies. For example, N-body simulations by 
Nipoti et al. (2006) revealed that disipationless collapses 
in preexisting dark matter haloes would naturally pro- 
duce galaxies which resemble ellipticals galaxies with de- 
pleted cores. These galaxies may also exist at z = 1.5 
to 2, however, it is unclear how these cores are protected 
against infalling high-density satellites in the absence of 
a central massive black hole. While the bulges of to- 
day's SOs do not all possess partially depleted cores, this 
does not exclude today's Sersic bulges from having ex- 
isted at z = 1.5 to 2. Indeed, as shown by MacArthur et 
al. (2009, and the references in Graham 2013) the bulk 
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Fig. .5. — Radial colour profiles for four of the five lenticular 
galaxies listed in Table [2] Colour information is not available for 
NGC 507. The inner R < 1" profile of NGC 3607 is excluded due 
to possible dust contamination. For NGC 6849, the B — R color 
profile is derived using the published surface brightness profiles 
presented in Reda et al. (2004). 

of most bulges' stars arc sufficiently old to have existed 
then. 

Using numerical simulations, Gocrdt et al. (2010) 
showed that the energy transferred from sinking massive 
objects can generate cores that are up to 3 kpc in sizes. 
Similarly, Martizzi et al. (2012a, see also Martizzi et al. 
2012b) used 500 pc resolution simulations and concluded 
that the combined effects of AGN feedback and inspiral- 
ing massive black holes create cores in luminous elliptical 
galaxies that are up to 8 - 10 kpc in sizes. However, these 
overly large cores reported by Gocrdt et al. (2010) and 
Martizzi ct al. (2012a) are in general incompatible with 
the majority of < 0.5 kpc cores observed in real galaxies 
(e.g., Trujil'lo et al. 2004; Richings et al. 2011; Dullo & 
Graham 2012). 

One can envisage that a core in a core-Sersic SO galaxy 
may be created through the 'dry' intermediate-mass 
merger (with mass ratio 1:4 to 1:7) of an SO galaxy 
with a smaller elliptical galaxy, such that the disk is 
not destroyed but the bulge grows while the black hole 
binary created during the merger ejects stars from the 
inner regions of the bulge. Eliche-Moral et al. (2012; 
see also Bekki 1998; Bournaud et al. 2005) has recently 
explored this merger scenario and concluded that 'dry' 
intermediate-mass mergers can give rise to remnants that 
are both photometrically and kinematically compatible 
with the observed SOs. It should, however, be noted that 
these studies did not explore the actual nuclear struc- 
ture of the merger remnants. Another possibility may 
be minor disk galaxy mergers which result in a rotating, 
disky merger remnant (Burkert & Naab 2005; Jesseit et 
al. 2005). Core-Sersic SO galaxies might also or instead 
be assembled via an inside-out mechanism: an early ma- 
jor merger might create an elliptical galaxy with a de- 
pleted core, around which a disk is subsequently built up 
over time via 'dry' minor mergers and/or smooth cold 
gas accretion. The assumption here is that the accretion 
process is very gentle and slow, enabling the cold gas 
to assemble into a star- forming disk (e.g., Steinmentz & 
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TABLE 3 
Core-Sersic galaxy data. 



Galaxy 


M/Lv 


log [AU/Mq) 


Be 1 


og {Laef/LQ,v) 


log (Maef/Me) 


log (Mbh/Mq) 








(kpc) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


NOO 0507 


5.5 


11.70 


1.65 


8.34 


9.08 


9.22 ± 0.39 


NGC 2300 


5.0 


11.16 


1.02 


7.81 


8.51 


8.86 ± 0.39 


NGC 3607 


5.7 


11.31 


0.84 








NGC 3706 


5.2 


11.60 


9.18 


8.30 


9.01 


8.93 ± 0.39 


NGC 4382 


2.6 


10.90 


0.99 


7.68 


8.09 


7.98 ± 0.38 


NGC 6849 


2.4 


11.32 


2.98 


7.97 


8.35 


8.04 ± 0.43 



Notes. — Col. (1) Galaxy name. Col. (2) V^-band stellar mass-to-light (M/L) ratio. Col. (3) Stellar 
mass of the bulge obtained using the bulge magnitude (Table [T} and the stellar mass-to-light (M/L) 
ratio (col. 2) (for the elliptical galaxy NGC 3706 the contribution of the galaxy's nuclear stellar ring has 
not been subtracted). Col. (4) Major-axis half-light radius of the bulge. Col. (5) Central luminosity 
deficit in terms of V-band solar luminosity. Col. (6) Central stellar mass deficit obtained using col. 
(2) and col. (5). Col. (7) SMBH mass predicted from the M — a relation presented in Graham et al. 
(2011). We use Eq. 4 from Graham et al. (2011) and the c- values in Table [T] and a 10% uncertainty on 
tr to estimate the error on SMBH mass. We adopt a barred morphology to estimate the mass of the 
black hole in NGC 6849. 



Navarro 2002; Birnboim & Dekel 2003; Keres et al. 2005) 
and also preventing it from falling towards the centers of 
the bulges, which otherwise would fuel a central burst of 
star formation and replenish the depleted cores. 

These later core-Sersic galaxy formation scenarios, 
each involving merger events, may be in agreement with 
the properties of our high luminosity SO bulges. This is in 
contrast to models that produce SOs by fading the disks 
of spiral galaxies through processes such as ram pressure 
stripping (e.g., Gunn & Gott 1972) — which is not to say 
that such processing does not additionally occur. 

6.2. Bulges of today's SOs versus compact high redshift 

galaxies 

Recent observations have revealed that the sizes of 
compact massive galaxies at redshift of z '--^ 1.5 to 2 (hav- 
ing effective radii i?e < 2 kpc) are up to a factor of ~ 5 
smaller than the present day elliptical galaxies of compa- 
rable stellar mass (Daddi et al. 2005; Trujillo et al. 2006; 
Damjanov et al. 2009; Saracco et al. 2011; Szomoru et 
al. 2012, and references therein). While today's ellipti- 
cal galaxies are widely regarded as the descendants of 
these compact high redshift galaxies (e.g., Hopkins et 
al. 2009), Graham (2013, his Fig. 1) showed that mas- 
sive local disk galaxies can also be alternative candidates 
given the high stellar density and compactness of their 
bulges. Comparing the properties of our massive present 
day bulges with compact high redshift galaxies, we plot 
the size-mass (Fig. [5^) and size-(Scrsic index) (Fig. ^jp) 
diagrams for a sample of 106 galaxies. 101 of these ob- 
jects are compact massive quiescent galaxies at redshift 
z = 0.2 to 2.7 taken from Damjanov et al. (2011, selected 
from their Table 2), while five are the bulges from local 
SOs (Table [2] and [S]). This figure clearly shows that the 
location of the bulges of our five core-Sersic SOs coincides 
with that of the high-redshift compact galaxies in both 
the Re — A'/h. and Re — n diagrams. 

This overlap suggests that the compact high redshift 
galaxies are the progenitors of (the bulges of) massive 
present day SOs (e.g.. Button et al. 2012; Graham et 
al. 2013). Indeed, the inside-out growth mechanism de- 
scribed at the end of Section 16. 1[ which could drive the 
build up of the core-Sersic SOs (Section [6.3p . would also 
eliminate the difficulties encountered in trying to explain 
the size evolution of these compact high redshift galaxies 



into much larger, modern day elliptical galaxies. Large 
elliptical galaxies already exist at z w 1.5 to 2 (e.g., Bruce 
et al. 2012). We further note that van der Wei. (2011, 
see also Chevance et al. 2012) wrote that the majority 
of compact high redshift galaxies have small undevel- 
oped disks. In addition, Poggianti et al. (2012) found 
that 4.4% of their total low redshift galaxy population 
had sizes and mass densities comparable to the compact, 
massive high redshift galaxies, with 70% of these com- 
pact low redshift galaxies found to be SOs. Given the 
results in Graham (2013, his Fig. 1), a comparison of lo- 
cal SO galaxy "bulge" sizes and densities, rather than the 
entire "galaxy" sizes and densities, should yield a higher 
percentage match. It seems plausible that minor merg- 
ers and cold gas accretion (e.g., Conselice et al. 2012), 
in a preferred plane due to known cosmological stream- 
ing/feeding paths, may transform some of the compact 
high redshift galaxies into modern SOs by the creation of 
a younger surrounding disk. This process eliminates the 
difficulties that arise when trying to evolve the compact 
high-z galaxies into local elliptical galaxies, such as the 
shortage of satellites required to produce the necessary 
expansion via many minor mergers (Trujillo et al. 2012). 
It is also consistent with cosmological models (e.g., Stein- 
metz & Navarro 2002), and explains many properties of 
early-type galaxies. 

6.3. The role of galaxy environment 

Several studies have highlighted that environment 
plays a major role in the formation of SO galaxies (e.g.. 
Dressier 1980; Dressier et al. 1997; Couch et al. 1998; 
Fasano et al. 2000; Poggianti et al. 2009; Wilman et al. 
2009; Just et al. 2010; Bekki & Couch 2011). Three 
of our five SO galaxies (NGC 507, NGC 2300 and NGC 
3607) reside in X-ray bright galaxy groups, while NGC 
6849 is an isolated galaxy, and NGC 4382 is a member 
of the Virgo cluster. 

NGC 6849, the isolated core-Sersic lenticular galaxy, 
may have had its bulge built through an early violent 
'dry' major merger, while subsequent late accretion of 
gas and stars built up its disk — the picture in Graham 
(2013), see also Conselice et al. (2012). This hierarchi- 
cal inside-out growth scenario is supported by its radial 
colour profile plotted in Fig. [51 The B — R colour gradi- 
ent reveals that the galaxy becomes progressively bluer 
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Fig. 6. — Major-axis effective radius Re plotted against (a) stellar mass M*, and (b) Scrsic index n for a sample of 106 galaxies. 101 
of these objects are massive compact galaxies, in the redshift range 0.2 < z < 2.7, published by Damjanov et al. (2011, their Table 2 
excluding galaxies (i) with no reported n values, (ii) fitted with the de Vaucoulers {R^^'^) model, and (iii) modelled by setting n to some 
constant value). The remaining five are local SO bulges from this study ( Table [2l and [3ll . Median uncertainty from Damjanov et al. (2011) 
for the 101 galaxies is shown by the error bar at the top of each panel. 



towards larger radii (> 10") where the disk dominates. 
Arnold et al. (2011) proposed such a two phase assembly 
mechanism to explain the isolated field galaxy NGC 3115 
(except for its dwarf companion) based on their analysis 
of its kinematics and metallicity. Rcda et al. (2004) also 
found that mergers are a dominant formation path to 
forming some isolated early- type galaxies. In passing we 
note that van den Bergh (2009) alternatively suggested 
that disk gas ejection by AGN might be a dominant pro- 
cess that transforms isolated spirals galaxies into SOs. 
While by itself this process docs not appear capable of 
generating a depleted stellar core in galaxies, it may have 
transformed some spiral galaxies into lenticular galaxies. 

The 'dry' major merger scenario for the bulges of bright 
lenticular galaxies in groups may be naturally consistent 
with the observed morphology density relation for the 
galaxy groups. Studies have shown that galaxy-galaxy 
interaction and merging play an important role in trans- 
forming group spirals into SOs and giving the observed 
SO fractions in groups (e.g., Postman & Gellerl984; Hels- 
don & Ponman 2003; Just et al. 2010; Bckki & Couch 
2011). Our three X-ray bright group galaxies are the 
brightest central galaxies from their respective groups. 
As pointed by Helsdon & Ponman (2003 and references 
therein) these central SO galaxies are assumed to be the 
byproducts of an earlier merger activity in the collaps- 
ing group core plus some recent merger and accretion 
events. Interestingly, as shown in Fig. [Sj the outskirts of 
these galaxies (colour information was not available for 
NGC 507) are relatively bluer than their centers. This 
result further strengthens the view that bulges in lumi- 
nous group SOs grow by mergers while their disks are 
gradually built and enhanced via latter gas accretion. 

For NGC 4382, our fifth lenticular galaxy and the only 
cluster galaxy in our sample, our analysis of its bright- 
ness profile and thus the mass deficit (Figs. |3] and 2]) sug- 
gests a 'dry' merger event as a formation path. In con- 
trast, processes such as ram pressure stripping (Gunn & 



Gott 1972; Quilis et al. 2000), galaxy harassment (Moore 
et al. 1996) and strangulation (Larson 1980; Bekki et 
al. 2002) are often noted as plausible formation mecha- 
nisms for cluster lenticular galaxies, rather than galaxy 
mergers, due to the high cluster velocity dispersion and 
high intracluster medium density. We do not deny these 
mechanisms, only that they alone cannot account for par- 
tially depleted cores. For example, Vollmer et al. (2008a, 
2008b) present wonderful observational evidence of the 
ram pressure stripped Virgo spiral galaxies NGC 4501 
and NGC 4522, while Merluzzi et al. (2012) describe 
in detail the stripping process of a galaxy in Abcll 3558. 
The merger event in NGC 4382's past may have occurred 
prior to it entering the cluster. Coupling N-body simu- 
lations with a semi-analytic formation model, Okamoto 
& Nagashima (2001) also reported that the fractions of 
cluster SO galaxies produced by major mergers are signif- 
icantly smaller than the observed fractions. In contrast, 
Bekki (1998) proposed unequal galaxy merging between 
two spirals as a dominant formation origin for cluster 
lenticular galaxies (see also Burkcrt & Naab 2005 and 
Jesseit et al. 2005). NGC 4382 is situated in the out- 
skirts of the Virgo cluster, a location where mergers are 
likely to occur, and it displays stellar shells in its im- 
age (Schweizer & Seitzer 1992). These appear to suggest 
that NGC 4382 might have a merger related origin (al- 
though see Chung et al. 2009). The g - z = 1.95 bulge 
colour of NGC 4382 (Fig. [5]) is however bluer than the 
typical g — z = 1.56 (AB mag) colour quoted for ellip- 
ticals (Fukugita et al. 1995). If the core in NGC 4382 
was formed from the inspiral of supermassive black holes 
in a relatively gas-free merger event, then the progenitor 
stars which make the bulge were not old, suggestive of 
ongoing core formation in some lenticular galaxies rather 
than all being formed at rcdshifts beyond 1.5 to 2 prior 
to the subsequent accretion of a disk. 

7. CONCLUSIONS 
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We have used the IRAF ellipse task to derive the 
major-axis surface brightness profiles and isophotal 
parameters for six early-type galaxies observed with 
the high-resolution HST WFPC2 and ACS cameras. 
While one of these turned out to be an elliptical galaxy, 
which we modelled with a core-Sersic profile, we have 
modeled the surface brightness profiles of the remaining 
five lenticular galaxies using a core-Sersic model for the 
bulge plus an exponential model for the disk. This is 
the first time this has been done for disk galaxies with 
depleted cores. In our analysis, we have additionally 
accounted for the bar component in one of these disk 
galaxies by using the Ferrers (1877) function. Our 
primary conclusions are as follows: 

1. The core-Sersic bulge plus exponential disk model 
gives an accurate description to core-Sersic lenticular 
galaxy light profiles. The rms residual scatter of the fits 
are < 0.05 mag arcsec"^. 

2. The Sersic index n is ^ 3 for the bulges of our core- 
Sersic lenticular galaxies, whereas n > 3 for the brighter 
core-Sersic elliptical galaxies in Dullo & Graham (2012). 

3. The core 'break radii' range from 24 to 102 pc. 

4. We have measured central stellar mass deficits in 
four of the luminous "core-Sersic" lenticular galaxies, 
finding Mdof ^ 0.5 — 2 Mbh, in agreement with previ- 
ous core-Sersic analysis of elliptical galaxies. (NGC 3607 
was excluded because its dusty nuclear spiral compro- 
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mises the recovery of its core parameters.) 

5. Our results tentatively suggest that, regardless of 
their environments, core-Sersic lenticular galaxies could 
be assembled in two stages: an earlier violent 'dry' ma- 
jor merger process involving massive black holes, which 
forms the bulge component, followed by subsequent disk 
formation through minor mergers and/or very gentle cold 
gas accretion about a preferred plane. Dry intermediate- 
mass ratio mergers of SOs with smaller BH-hosting ellip- 
ticals could also account for the buildup of some core- 
Sersic SOs. 

6. The location of the bulges of our five SO galaxies, 
including the SO NGC 3607 possibly with out a depleted 
core, in the mass-size and mass- (Sersic index) diagram 
coincides with that of the compact early-type galaxies 
seen at redshift of z ~ 0.2 to 2.7. This could be an indi- 
cation that today's massive bulges may be descendants 
of the compact high redshift early-type galaxies. 
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9. APPENDIX A 

9.1. Notes on individual galaxies 

In this section we review several relevant features associated with each galaxy in our sample. We comment on the 
X-ray properties as this will be relevant to our later discussion about formation scenarios and environment. 

9.1.1. NGC 507 

NGC 507 is the brightest galaxy in the nearby, poor NGC 507 group; it is also one of the brightest known X-ray 
early-type galaxies. Laurikainen et al. (2010) found evidence of a weak bar in their Kg-hand image. However, they 
fitted a Sersic bulge plus an exponential disk to their deep sub-arcsec resolution Ks-hand data because the bar was 
too weak to be included in the fit. Our high-resolution HST/F555W surface brightness profile for this galaxy is well 
fit by the core-Scrsic bulge -I- exponential disk model with an rms residual of 0.027 mag arcsec"^ (Fig. [3)). 

9.1.2. NGC 2300 

NGC 2300 belongs to the poor NGC 2300 group. Mulchaey et al. (1993) reported the presence of hot diffuse 
intragroup gas close to NGC 2300 using ROSAT. Although several previous studies (e.g., Huchtmeier 1994, Sandage 
& Bedke 1994) have classified this galaxy as an elliptical, as shown in Fig. [3l the luminosity profile is well fit by 
the core-Sersic bulge -|- exponential disk model (see also Laurikainen et al. 2010 for their Sersic -I- exponential fit to 
ground-based data that does not resolve the depleted core.) 

9.1.3. NGC 3607 

NGC 3607 is an X-ray luminous (e.g., Terashima et al. 2000), bright lenticular galaxy (e.g., de Vaucouleurs et al. 
1991; Afanasiev & Silchenko 2007; Laurikainen et al. 2010; Gutierrez et al. 2011; although it is classified as an elliptical 
in Giiltekin et al. 2009) located in the Leo II group. It has a dust obscured central region as can be seen from the 
residual image in Appendix B, Fig. HI We constructed an image model using the IRAF task bmodel and subtracted 
it from the raw F-band image to obtain this residual image. The effect of this dusty nuclear spiral can also be seen in 
the position angle, ellipticity and isophote shape parameter (B^) profiles (Fig. [5]). In addition to using this to guide 
the careful dust masking procedure, performed prior to running ellipse to obtain the y-band light profile beyond 1", 
we have analysed the galaxy's near-infrared HST/NICMOS F160W (_H-band) image which is less affected by dust. As 
noted in Section 2, the inner {R < 1") V-band profile of this galaxy, taken from Lauer et al. (2005), is derived from a 
PSF-deconvolved image. Therefore, the apparent 1".31 core of the F-band profile fit is most likely an artifact due to 
interference from the dusty nuclear spiral disk with the deconvolution routine (as noted already by DuUo & Graham 
2012, and references therein). 

9.1.4. NGC 3706 

In agreement with the ground-based work by Laurikainen et al. (2010), the analysis of NGC 3706's HST/F555W 
light profile suggests that it is likely an elliptical galaxy rather than a lenticular galaxy. The unsharp-masked image 
(Appendix B, Fig. [8]), the position angle twist, the disky and very flat ('pointy') nature of the isophotes from the 
ELLIPSE fit in the region i?= 0".2 — 1".0 (Fig. [2]) reveal the presence of an edge-on nuclear ring of stars in this galaxy 
(Lauer et al. 2002; Kandrup et al. 2003). We model the host galaxy light profile after excluding the region affected by 
this additional ring of star light (Fig. [3]). 

9.1.5. NGC 4382 

NGC 4382 is a peculiar, fast rotating galaxy (Emsellem et al. 2007) in the Virgo cluster which is classified as an SO 
in the RC3. Based on the photometry and brightness profile analysis, but not the rotational evidence, Kormendy et 
al. (2009) reported an absence of a large-scale stellar disk in this galaxy, and identified the extra light at large radii 
to be a feature associated with a recent merger (although see Chung et al. 2009), hence, they classified the galaxy as 
E2. In contrast, Laurikainen et al. (2011) have classified the galaxy as an SO based on the detection of dispersed spiral 
arm segments, although they remained uncertain about the exponential nature of the disk light distribution. In Fig. [31 
we show that the galaxy light profile is indeed well fit by the core-Sersic bulge -I- exponential disk model with an rms 
residual of only 0.016 mag arcsec"^. This galaxy has positive B4 values, i.e. disky isophotes where R > 1". We also 
note that we measure the galaxy's bulge-to-disk (B/D) flux ratio to be ^0.67 (B/T=0.40), revealing the existence of 
an appreciable fraction of the galaxy light in the large-scale disk. Dullo & Graham (2012, their Fig. 7) show how this 
large-scale disk results in a different "galaxy" Sersic index, from the single-component fit in Ferrarese et al. (2006), 
compared to the "bulge" Sersic index. Not fitting for the outer disk results in a 'galaxy' Sersic index which is larger 
than the "bulge" Sersic index and thus over-predicts the central mass deficit. 



Depleted cores in lenticular galaxies 



15 




Fig. 7. — Double Sersic model fits to the surface brightness profile of NGC 1374 (inner exponential + Sersic), NGC 4473 (inner 
exponential + Sersic) and NGC 5576 (inner Sersic + outer Sersic), see the text for details. The ~ 100" geometric mean light profiles for 
NGC 1374 and NGC 4473 were taken from Turner et al. (2012) and Ferrarese et al. (2006), respectively. For NGC 5576, we use the Lauer 
(2012) composite ( deconvolved HST data of Lauer et al. 2005 at /? < 10" + ground-based data at ij > 10") major-axis light profile 
which extends out to ~ 100". In contrast to Lauer (2012, his Fig. 1) and Lauer et al. (2005), our modelling (as in DuUo & Graham 2012) 
reveals the absence of depleted cores in these galaxies. The rms residuals. A, are shown in each panel. 

9.1.6. NGC 6849 

NGC 6849 is classified as an isolated, barred lenticular (SBO) galaxy in the RC3. It is also classified as an elliptical 
galaxy in the ground-based studies of Saraiva et al. (1999) and de Souza, Gadotti & dos Anjos (2004). In contrast 
to the findings of Saraiva et al. (1999), the F555W/775'T' unsharp-masking (Appendix B, Fig. [8|) appears to suggest 
the presence of a bar in this galaxy. Fig. [5] illustrates that the ellipticity of this galaxy increases towards larger radii, 
which makes it the flattest in our sample. As shown in Fig. [3l the 3-component bulge-bar-disk luminosity profile is 
well fit by the the core-Sersic bulge plus Ferrers bar plus exponential disk decomposition model with an rms residual 
of ^0.04 mag arcsec^^. 

9.2. Comments on Dullo & Graham (2012) 

Commenting on Dullo & Graham (2012), Lauer (2012) presented Sersic model fits to the outer regions of composite 
light profiles (using the inner 10" HST profiles from Lauer et al. 2005 combined with ground-based profiles out to 
^100" taken from the literature) for three galaxies which Dullo & Graham (2012), and others, had shown not to 
possess depleted cores: NGC 1374, NGC 4473 and NGC 5576. Lauer's (2012) fits suggested the presence of depleted 
cores relative to his outer Sersic model. He therefore claimed that Dullo & Graham (2012) had misclassified these 
three galaxies as coreless Sersic galaxies due to an incorrect application of the Sersic model to radially limited 10" light 
profiles taken from Lauer et al. (2005). However, there were already ~100" high-resolution HST profiles published 
for these galaxies that had been well fit using Sersic models with no central light deficit (Turner et al. 2012, NGC 
1374; Ferrearese et al. 2006 and Kormendy et al. 2009, NGC 4473; Trujillo et al. 2004, NGC 5576). Fig. [7] shows 
that the light profiles for these galaxies, covering a large radial extent, are described by the sum of two Sersic profiles 
without any partially depicted core. The reason for these two-component models is detailed below. Basically, taking 
into account additional information, such as the presence of tidal material and/or kinematic substructure can help to 
identify the required components of a fit. In general, we recommend showing the residuals about one's fitted model, 
as they can reveal if the large-scale curvature in the radial light distribution is well matched, and we advise against 
subjectively restricting the fitted radial range. 



9.2.1. NGC 1374 

Turner et al. (2012) fit a PSF-convolved Sersic model to this galaxy, with no evidence for a depleted core. In Dullo 
& Graham (2012) we revealed that there is actually an additional component at the center of this galaxy. In Fig. [7^, 
we show our (PSF convolved) nuclear disk -I- Sersic model fit to the light profile of NGC 1374 taken from Turner et al. 
(2012), see also the fit presented in Dullo & Graham (2012) to the inner ~ 10". Some small (in amplitude), large-scale 
residual is present in Fig. [7^. The rotation curve in Graham et al. (1998, see also Longo et al. 1994) revealed that 
this galaxy rotates out to 40", suggesting the presence of an additional component which may be the cause of some of 
the residual structure about our exponential (for the nuclear disk) plus n = 4.3 Sersic (for the underlying host galaxy 
light) model. 
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9.2.2. NGC 4473 

Ferrarcse ct al. (2006) fit a PSF-convolvcd Scrsic model to NGC 4473, with no depleted core. Kormendy et al. (2009) 
identified a central light excess in this galaxy, above their adopted Scrsic {n « 4) fit, and they remarked that the extra 
light is due to a known counter-rotating embedded stellar disk (Emsellem et al. 2004; Cappellari & McDermid 2005; 
Cappellari et al. 2007). Fig. [7}d shows our convolved exponential disk + Sersic bulge model for this galaxy's light 
profile taken from Ferrarese et al. (2006). We fit an exponential n = 1 function to the counter rotating disk, and find 
that a Sersic n ^ 3 model fits the underlying galaxy light. Furthermore, Pinkney et al. (2003) remarked that this 
galaxy has unusual parameters from their Nuker model fit. They found that it has i) the smallest a value of their 
sample (i.e. a broad transition region between the inner and outer Nuker model power-law profiles) and ii) a very 
steep /3 value (the steepest from all the galaxies modelled by Byun et al. 1996). As warned in Graham et al. (2003), 
this is what one expects when fitting the Nuker model to what is actually a Sersic profile with a low value of n and 
no depleted core. Pinkney et al. (2003) additionally remarked that the absolute magnitude of this galaxy is consistent 
with the "power-law" galaxies, i.e. those without depleted cores. 

9.2.3. NGC 5576 

Trujillo et al. (2004) modelled the deconvolved 100" light profile of NGC 5576 with the Sersic model, finding no 
depleted core. Kim et al. (2012, their Fig. 5) showed that NGC 5576 and the barred SO NGC 5574 are an interacting 
pair (see also Tal et al. 2009) having a long tidal tail, which extends out to 60 kpc in their 3.6 /.m image. They showed 
the presence of tidal disturbances in the outskirts (from 1' to 2') of NGC 5576 because of this interaction. In Fig. [7J:, 
we therefore fit a double Sersic model to the light profile of NGC 5576 presented in Lauer (2012)0. We find that the 
central galaxy light distribution is well fit with an inner n ~ 3.5 Sersic model, and the extra light, which is likely to be 
due to the re-distribution of material from the minor to the major-axis as a result of the tidal interaction, is described 
by an outer n '-^ 1.2 Sersic model. Kim et al. (2012) also fit a double Sersic model (an inner n = 3.45 Sersic model 
plus an outer n = 1.48 Sersic model) to their light profile which is sampled from R = 11". 8 to 247". 5. In general 
agreement with this and the inner n ~ 4.5 Sersic model fit presented in Trujillo et al. (2004), our n ~ 3.5 Sersic fit to 
this galaxy's inner light distribution (Fig. [T]:) not only shows the absence of a central luminosity deficit in this galaxy, 
but is also consistent with its magnitude My = —21.11 mag and velocity dispersion a = 171 km s^^ (see Graham et 
al. 2001, their Fig. 13; DuUo & Graham 2012, their Fig. 14). This is in contrast to the single-component n '--^ 8 Sersic 
fit by Lauer (2012, then Fig 1). 

We also note that, contrary to the claim by Lauer (2012), DuUo & Graham (2012, their section 4) explicitly noted 
that the spurious downward departure in the inner light profile of NGC 4486B (relative to the inward extrapolation of 
its outer Sersic profile) is due to the presence of a double optical nucleus, rather than a typical depleted core. Coupled 
with NGC 4458, NGC 4478 and NGC 7213, we therefore have that seven of the 39 galaxies (i.e. -18%) classified 
as "core" galaxies according to the Nuker model analysis by Lauer et al. (2005), do not have central stellar deficits 
relative to their outer Scrsic profiles and are therefore not identified as "core-Scrsic" galaxies. Such 'mismatches' will 
be most prevalent among galaxy samples fainter than AIb ~ —20.5 mag that contain spheroids with low Scrsic indices 
and thus have relatively flat inner cores that are not depleted of stars. 

The ATLAS 3D galaxy sample used by Lauer (2012) is dominated by early-type galaxies with magnitudes up to 
~3 mag fainter than Mb = — 21 mag — many of which are two-component lenticular galaxies (Emsellem et al. 2011). 
As such, application of the Nuker model may identify many 'cores' (7 < 0.3) among these galaxies which do not 
actually possess partially-depleted cores relative to their spheroid's outer Sersic profile. The Nuker model classification 
may therefore substantially blur the actual connection between the dry merger scenario, as described by Faber et al. 
(1997), and the existence of depleted cores from coalescing black holes. We thus advocate the core-Scrsic model and 
point readers to additional important reasons for this that are discussed in Graham et al. (2003) and DuUo & Graham 
(2012). 

10. APPENDIX B 

Fig. [8] shows distinct features in NGC 3607 (left), NGC 3706 (middle) and NGC 6849 (right). 



® Lauer (2012) noted that his ground-based profile for NGC 5576 , i n ■ 1 j j u tu ti, 

■jju ni-i JO iM i 1 /-.r,r.o\ 1 ■ 1 was uot actuallv mcludcd by those authors, 

was provided by Michard & Marcheal (1993), however this galaxy •' 
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Fig. 8. — Left: Residual image of NGC 3607 obtained by subtracting our symmetrical image model from the original V^band image of 
the galaxy observed using the HST/WFPC2 camera. A dust disk is seen in the central region of the galaxy, see the Appendix A text for a 
further description. Middle: H5T WFPC2/F555W PC 15"xl5" unsharp masked image of NGC 3706 showing an edge-on nuclear stellar 
ring (sec Fig. [3] and Laucr ct al. 2002, their Figure 1). Right: HST WFPC2/F555W PC 15"xl5" unsharp masked image of NGC 6849 
showing an inner bar. For NGC 3607, north is in the direction of the arrow. North is up for NGC 3706 and NGC 6849. 



